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Integrable  Models  of  Shallow  Water  Waves 

by 

Harvey  Segur*,  Allan  Finkel**, 
and  Hilda  Philander* 


Abstract: 

The  Korteweg-de  Vries  and  the  Kadomtsev-Petviashvil i  (KP)  equations  both 
model  the  evolution  of  relatively  long  water  waves  of  moderate  amplitude, as 
they  propagate  in  shallow  water.  Both  equations  are  completely  integrable. 
In  this  paper  we  review  the  derivation  of  each  equation  as  an  approximate 
model  of  shallow  water  waves,  and  compare  their  solutions  with  some  of  the 
experimental  observations  of  waves  in  shallow  water.  We  also  describe  in 
detail  the  family  of  doubly  periodic  KP  solutions.  These  are  the  natural 
two-dimensional  generalizations  of  cnoidal  waves  in  one-dimension;  one  may 
think  of  them  as  describing  "typical"  patterns  of  nonlinear,  two-dimensional 
waves  in  shallow  water. 
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08540 
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I.  Introduction 


Two  equations  that  have  been  studied  intensively  in  recent  years  are  the 
Korteweg-de  Vries  (KdV;  1895)  equation. 


u  +  6uu  +  u  =0. 

t  x  XXX  » 


(1) 


and  a  generalization  of  it  due  to  Kadomtsev  and  Petviashvili  (KP;  1970), 


(u  +  6uu  +  u  )  +  3u  =0. 

v  t  x  xxx'x  nn 


(2) 


Most  of  this  interest  is  due  to  the  remarkable  fact  that  each  equation  can  be 
solved  exactly  as  an  initial-value  problem  by  a  method  now  known  as  the 
Inverse  Scattering  Transform.  This  method  was  first  discovered  for  the  KdV 
equation  on  -«  <  x  <  «,  in  the  famous  papers  of  Gardner,  Greene,  KrusKal  and 
Miura  (1967,  1974).  The  corresponding  work  for  (1)  on  a  periodic  interval  was 
published  by  several  people  during  1974-1976  [Novikov  (1974),  Dubrovin  and 
Novikov  (1974),  Dubrovin  (1975),  Lax  (1975),  Its  and  Matveev  (1975),  McKean 
and  van  Moerbeke  (1975),  McKean  and  Trubowitz  (1976),  Dubrovin,  Matveev  and 
Novikov  (1976)].  For  the  KP  equation  on  -®  <  x»  h  <  *,  a  method  of  solution 
was  given  very  recently  by  Ablowitz,  Bar  Yaacov  and  Fokas  (1982;  cf.  Ablowitz 
and  Fokas,  these  Proceedings). 

As  we  shall  see  below,  it  happens  that  both  (1)  and  (2)  also  model  the 
evolution  of  water  waves  of  moderate  amplitude  as  they  propagate  in  one 
direction  in  relatively  shallow  water.  In  physical  terms,  the  KdV  equation 
arises  if  the  waves  are  strictly  one-dimensional  (i.e.,one  spatial  dimension 
plus  time),  while  the  KP  equation  arises  if  they  are  only  nearly 
one-dimensional.  Because  both  equations  are  completely  integrable,  they 
provide  very  precise  predictions  about  the  evolution  of  water  waves  under 
appropriate  conditions. 

This  paper  has  two  objectives.  The  first  is  to  examine  the  validity  of 
(1)  and  (2)  as  models  of  water  waves,  by  comparing  their  solutions  with  some 
of  the  available  experimental  data.  The  comparison  given  here  will  be  brief. 
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but  much  more  detailed  verifications  of  (1)  have  been  made  elsewhere  (e.g., 
Hammack  and  Segur,  1974,  1978).  The  second  objective  is  to  describe  in  detail 
a  special  family  of  solutions  of  (2),  which  are  periodic  in  each  of  two 
independent  variables.  It  appears  that  these  doubly-periodic  solutions  may 
have  great  practical  importance  as  models  of  long  water  waves.  Among  other 
things,  they  describe:  (i)  the  nonlinear  interaction  of  two  trains  of  finite 
amplitude  waves  in  shallow  water;  (ii)  the  reflection  of  a  train  of 
finite-amplitude  shallow-water  waves  from  a  vertical  wall  (by  replacing  an 
appropriate  line  of  symmetry  with  the  wall);  ( i i i )  the  reflection  of  such  a 
wave  train  by  a  change  in  bottom  topography;  or  (iv)  "typical" 
finite-amplitude,  short-crested  waves  in  shallow  water.  In  this  last  sense, 
they  are  the  natural  generalizations  to  two  dimensions  of  cnoidal  waves  in  one 
dimension. 

II.  Derivation  of  the  Equations 

i 

The  classical  problem  of  water  waves  is  to  find  the  irrotational  motion 
of  an  inviscid,  incompressible,  homogeneous  fluid,  subject  to  a  constant 
gravitational  force  (g).  The  fluid  rests  on  a  horizontal  and  impermeable  bed 
of  infinite  extent  at  z  =  -h  and  has  a  free  surface  at  z  =  ?(x,y,t);  see 
Figure  1.  In  this  derivation  we  neglect  the  effects  of  surface  tension  at  the 
free  surface,  although  it  can  be  included  without  difficulty  (e.g.,  see  Ch.  4 
of  Ablowitz  and  Segur,  1981). 

INSERT  FIGURE  1  ABOUT  HERE 

The  fluid  has  a  velocity  potential,  41,  which  satisfies 

v2  4  =  0,  -h<z<c(x,y,t);  (3) 

(irrotational  motion  of  an  incompressible  fluid).  It  is  subject  to  boundary 
conditions  on  the  bottom,  z  =  -h: 

4Z  =  0,  (4) 
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(impermeable  bed);  and  along  the  free  surface,  z  =  5: 


5ts 

Dt  " 


tt  +  Mx  +  *ycy  =  *z 


(5) 


(kinematic  condition); 


*t  +  9C  + 


1 

2 


(6) 


(dynamic  condition). 

Boundary  conditions  in  (x,y)  and  initial  conditions  also  are  required. 
If  the  waves  in  question  are  isolated,  then  v$  and  5  should  vanish  as  (x  + 
y2)-M*>.  In  other  problems,  periodic  boundary  conditions  in  x  and  in  y  may  be 
relevant. 

This  problem,  first  posed  by  Stokes  (1847),  remains  unsolved.  Tomake 
further  progress,  we  impose  additional  assumptions  on  the  solutions  of 
(3 )- (6) .  The  first  such  assumption  is  that  the  wave  amplitudes  should  be 
small.  If  we  interpret  "small"  to  mean  infinitesimal,  then  we  may  linearize 
(3)-(6)  about  =  0,  c  =  0,  and  seek  solutions  of  the  linearized  equations 

proportional  to  exp  (i (kx+my-uit)}  (e.g.,  see  Lamb,  1932,  §§  228,  266,  267). 
The  result  is  the  linearized  dispersion  relation. 


u)  *  g<  tanh  xh 


(7) 


2  2  2 

where  x  =  k  +m  .  From  this  one  computes  the  group  velocity  and  shows  that 
the  linearized  problem  is  dispersive  at  most  wave  numbers,  but  not  as  x+O, 
(i.e.,  long  waves,  or  shallow  water  waves),  where  it  is  only  weakly 
dispersive.  Both  (1)  and  (2)  arise  as  models  of  the  water  wave  problem  in 
this  weakly  dispersive  limit,  xh«l. 
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To  derive  (1)  or  (2),  we  assume  that: 


(A)  wave  amplitudes  are  small. 


e  s  C  max/h  <<  1; 


(B)  the  relevant  length  scale  in  the  x-direction  is  much  longer  than  the 
fluid  depth  (i.e.,  shallow  water  waves), 

(kh)2  «  1; 

(C)  either  the  motion  is  strictly  one-dimensional  (for  KdV), 


m  =  0, 


(C')  or  it  is  nearly  one-dimensional  (for  KP) 


(?)’ 


«  l; 


(D)  All  of  these  effects  balance, 

KdV:  (kh)2  =  0(e),  KP:  (kh)2  =  0(e)  =  0(mh). 

These  assumptions  imply  a  certain  scaling  of  the  original  equations,  (see  (9) 
below),  which  may  then  be  solved  as  a  perturbation  series,  term-by-term  in  e. 
[We  omit  the  details  here,  but  the  entire  derivation  may  be  found  in  Chapter  4 
of  Ablowitz  and  Segur  (1981),  or  elsewhere.] 

At  leading  order,  the  equations  are  linear  (from  (A)),  nondispersi ve 
(from  (B) ) ,  and  one-dimensional  (from  (C)).  The  solution  is  simply  the 
solution  of  the  linear,  one-dimensional  wave  equation: 

t  (x.y.t)  =  eh  [f (x-/gh  t;y)  +  F(x+/gh  t;y )]  +  0(e2).  (8) 
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At  this  order,  every  wave  has  permanent  form,  not  because  it  is  a  soliton 
but  because  we  are  solving  the  one-dimensional,  linear  wave  equation. 

When  this  perturbation  expansion  is  carried  to  second  order,  the 
one-dimensional,  linear  wave  equation  has  homogeneous  (forcing)  terms 
representing  weak  nonlinearity,  weak  dispersion  and  weak  two-dimensionality. 
Each  of  these  effects  contributes  to  a  secular  term  at  second  order.  We  may 
eliminate  them  by  introducing  a  second  slow  time-scale  (T): 

r  =  /e  (x-/gh  t)/h,  1  =  /e  (x+/gh  t)/h, 


o  =  ey/h,  T  =  ^  /g/h  t. 


(9) 


and  requiring  that  the  right-running  waves  satisfy 


(2  fT  ♦  3ffp  +|frpr)r  +  fnn  =  0, 


'(10) 


which  may  be  rescaled  to  (2).  In  the  one-dimensional  case,  3^=0,  and  (10) 
leads  to  (1).  The  left-running  waves  must  satisfy  a  similar  equation. 


At  this  order,  no  secular  terms  arise  from  interactions  between  left-  and 
right-running  waves  provided  that  f  and  F  are  smooth,  and  that  JVdr,  JVdl 
remain  bounded.  For  periodic  initial  data,  the  latter  requirement  implies 
that 


(ID 


where  the  integrals  are  taken  over  one  period.  Because  the  original  equations 
are  gal ilean-invariant,  (11)  amounts  only  to  a  normalization. 


Some  observations  about  the  physical  meaning  of  (1)  and  (2)  may  be  made 
at  this  point. 
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(A)  The  time  in  (1)  or  (2)  has  the  physical  meaning  of  a  slow  time-scale. 
Concepts  like  solitons  that  are  implied  by  these  equations  have  physical 
meaning  only  on  this  long  time-scale. 

(B)  Our  original  equations,  (3)  -  (6),  are  themselves  only  approximately 
correct,  because  they  neglect  real  effects  like  viscosity.  Thus,  there 
is  a  short  (linear,  non-dispersi ve,  one-dimensional),  time-scale,  on 
which  all  waves  have  permanent  form,  from  (8).  This  is  followed  by  a 
longer  (KdV  or  KP)  time-scale  on  which  soliton  interactions  are  relevant. 
This  is  then  followed  by  an  even  longer  time-scale  on  which  other  effects 
(like  viscosity)  become  important. 

(C)  The  KdV  and  KP  equations  are  fully  nonlinear,  but  they  model  water  waves 
only  when  the  water  waves  are  weakly  nonlinear.  Fully  nonlinear  water 
waves  are  known  to  break,  as  solutions  of  (1)  or  (2)  do  not.  Similarly, 
(2)  models  water  waves  that  are  only  weakly  two-dimensional,  although 
this  restriction  is  not  evident  in  (2)  by  itself. 

(D)  Because  (1)  and  (2)  are  first-order  in  time,  neither  describes  the 
interactions  of  left-  with  right-running  waves.  However,  this  is  not 
because  such  interactions  were  not  admitted,  but  because  they  are  not 
important  on  the  time-scale  on  which  these  equations  apply.  Their 
interaction  is  given  by  (8),  to  leading  order. 

III.  Experimental  Evidence  of  Solitons 

The  consequences  of  the  theory  of  Inverse  Scattering  Transforms  for  the 
KdV  equation  are  as  follows.  (Details  may  be  found  in  Ablowitz  and 
Segur,  1981,  among  other  places.)  Let  uQ(x)  be  any  smooth  function  on  (-•,») 
that  vanishes  rapidly  along  with  its  derivatives  as  x+±“*  Then  (1)  has  a 
unique  solution  that  coincides  with  uQ(x)  at  t=o.  As  this  solution 

evolves  into  N  isolated  solitons,  ordered  by  amplitude,  followed  by  an 
oscillatory  wave-train  that  disperses  in  time  ("radiation").  When  the 
solitons  have  separated,  each  is  given  (locally)  by 


7 


(12) 


u(x,t)  =  2k2  sech2  {k(x-4k2T+x0)}. 

The  number  of  solitons  (N),  their  amplitudes  (kj,  j=l,  ....  N),  and  all  the 
other  details  of  the  long-time  solution  can  be  found  directly  from  u0(x). 
There  is  no  need  to  advance  (1)  numerically  in  time. 

Corresponding  to  (12)  is  a  water  wave,  whose  surface  elevation  is  given 
by 

+  0(e2).  (13) 

Similarly,  all  the  other  consequences  of  (1)  imply  predictions  about  real 
water  waves.  Hammack  and  Segur  (1974,  1978)  tested  several  aspects  of  this 
theory  by  comparing  it  with  laboratory  experiments.  Next  we  briefly  reiterate 
some  of  their  results,  to  give  the  reader  an  idea  of  the  validity  of  the  KdV 
equation  as  a  model  of  long  water  waves  of  moderate  amplitude.  (Experimental 
comparisons  of  (1)  with  water  waves  also  have  been  made  by  Zabusky  and 
Galvin  (1971),  and  by  Weidman  and  Maxworthy  (1978).) 

INSERT  FIGURE  2  HERE  OR  BELOW. 


/  2  \ 

c(x,t)  =  ek2h  sech2  < 

✓ek 

h 

x±/gh  ( 1  +  It  +  X0 

The  experiments  of  Hammack  and  Segur  were  conducted  in  a  wave  tank 
31.6  m.  long,  61  cm.  deep  and  39.4  cm.  wide.  As  shown  schematically  in  Figure 
2,  the  wave  generator  consisted  of  a  rectangular  piston  located  in  the  tank 
bed  adjacent  to  the  upstream  end  wall  of  the  tank.  The  piston  spanned  the 
tank  width,  and  was  61  cm.  long  for  the  experiments  we  will  discuss.  The 
time-history  of  its  vertical  displacement  was  prescribed  for  each  experiment. 

Wave  measurements  were  made  during  each  experiment  at  several  positions 
down  the  tank  using  parallel-wire  resistance  gauges.  In  the  experiments 
described  here  the  fluid  depth  (h)  was  5  cm.,  and  the  waves  were  measured  at 
x/h  =  0,  20,  180  and  400,  where  x  =  0  at  the  downstream  edge  of  the  piston. 
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Figure  3  shows  a  wave  generated  simply  by  raising  the  piston.  The  piston 
motion  was  fast  enough  that  the  shape  of  the  wave  at  x/h  =  0  is  effectively 
the  shape  of  the  piston;  (because  of  the  reflecting  wall  at  the  upstream  end 
of  the  piston,  the  wave  at  x  =  0  was  actually  twice  as  long  as  the  piston  and 
half  as  high  as  its  displacement). 

INSERT  FIGURE  3  SOMEWHERE  NEAR  HERE. 

On  a  short  time-scale,  according  to  (8),  this  wave  should  simply 
translate  with  speed  v^gF.  The  wave  measured  at  x/h  =  20  (Figure  3b)  fits  this 
description  approximately;  its  shape  is  basically  that  of  the  wave  at  x  =  0. 
(The  front  of  the  wave  is  to  the  left  in  these  figures,  and  a  wave  which 
translates  with  speed  /gF  shows  no  horizontal  displacement  in  succeeding 

frames). 

That  solitons  emerge  on  a  long  time  scale  may  be  seen  in  Figures  3c,  d. 
Solving  the  appropriate  scattering  problem  with  the  wave  measured  at  x  =0  as 
the  potential  yields  3  discrete  eigenvalues,  representing  3  solitons.  These 
correspond  to  the  3  positive,  more-or-less  permanent  waves  seen  at  x/h  =  180 
and  x/h  =  400.  According  to  (13),  these  waves  all  should  move  to  the  left  in 
these  figures,  since  their  speeds  all  exceed  »*gF.  That  they  do  not  is  a 
measure  of  the  effect  of  viscosity  in  these  experiments. 

Even  so,  we  assert  that  these  waves  are  solitons  on  the  basis  of  their 
shapes.  The  entire  profile  of  a  single  soliton  is  determined  from  (13)  once 
its  amplitude  is  known.  The  peak  amplitudes  of  the  first  two  waves  in 

Figure  3d  were  measured  and  the  dots  in  that  figure  represent  evaluations  of 
(13)  based  on  those  amplitudes.  The  agreement  with  the  measured  wave  shapes 
is  striking. 

The  results  shown  in  Figure  3  suggest  the  following  picture  of  long  water 
waves  of  moderate  amplitude. 

a)  There  is  a  short  (linear)  time-scale,  during  which  the  left-  and 

right-running  waves  separate  from  each  other. 
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b)  There  is  a  long  (KdV)  time-scale,  during  which  the  right-  (or  left-) 
running  waves  evolve  in  N  solitons  plus  radiation. 

c)  There  is  an  even  longer  viscous  time-scale,  during  which  the  energy  in 
these  solitons  is  gradually  dissipated.  Because  the  KdV  time-scale  is 
shorter,  however,  the  solitons  continually  readjust  their  shapes  and 
speeds  as  they  lose  energy  so  that  local ly,  as  in  Figure  3d,  they  look 
and  act  like  solitons. 

Several  other  experiments  were  performed  in  this  series,  to  test  other 
aspects  of  the  KdV  theory.  We  refer  the  reader  to  the  original  papers  or  to 
Chapter  4  of  Ablowitz  and  Segur  (1981)  for  more  details. 

Next  we  consider  the  KP  equation,  (2),  as  a  model  of  nearly 

one-dimensional  water  waves  of  moderate  amplitude,  propagating  in  shallow 

water.  Here  the  theory  is  still  incomplete,  and  we  are  aware  of  no  systematic 
experimental  study.  Consequently  we  are  forced  to  discuss  special  solutions 
of  (2),  and  fortuitous  experimental  observations. 

Obviously,  every  KdV  solution  also  solves  KP.  More  generally, 

Satsuma  (1976)  showed  that  (2)  admits  an  N-soliton  solution,  with  the  N 

solitons  traveling  in  N  different  directions.  The  two-soliton  formula  is 

u(x»n»x)  =  2  3^  *n  f,  (14) 


where 


f  =  1  +  exp  (4^)  +  exp  ($2)  +  exp  U^^+A), 


4> j  =  kj  (x  +  Pj1!  “  cjT)» 


,2  „  2 
=  ki  +  3Pi 


t 


exp  (A)  = 


(krk2)2  '  (prp2)2 

(kj+k2)2  -  <prp2)2 
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A  typical  solution  is  shown  in  Figure  4.  Far  from  the  interaction  region, 
each  wave  is  essentially  a  KdV  soliton,  but  traveling  at  an  angle  to  the 
X-axis.  The  interaction  of  the  two  waves  is  necessarily  nonlinear,  and  a 
phase  shift  of  each  wave  as  a  result  of  the  interaction  is  evident.  There  are 
two  phases  in  this  solution  and  two  spatial  coordinates,  so  unless  Pj=P2 
(14),  there  is  a  uniformly  translating  coordinate  system  in  which  this 
solution  is  stationary. 

FIGURES  4  &  5  GO  ON  THE  SAME  PAGE,  SOMEWHERE  NEAR  HERE 

How  well  does  this  two-sol iton  solution  of  (2)  predict  water  wave 
interactions?  No  quantitative  experimental  data  is  available  but  the 
photograph  in  Figure  5,  of  two  long-crested  waves  interacting  in  shallow  water 
off  a  beach  in  Oregon,  certainly  is  suggestive.  Each  of  the  two  waves 
apparently  is  part  of  a  train  of  periodic  waves  coming  in  from  deep  water,  but 
their  wavelengths  seem  to  be  long  enough  that  each  acts  like  a  solitary  wave 
in  shallow  water.  The  comparison  of  Figures  4  and  5  is  only  qualitative,  of 
course,  but  certainly  it  suggests  that  the  KP  equation  might  be  as  useful  for 
weakly  two-dimensional  water  waves  as  the  KdV  equation  is  for  one-dimensional 
waves. 

Because  (2)  only  models  weakly  two-dimensional  wave  interactions,  some 
oblique  interactions  are  not  predicted  by  it.  Miles  (1977a, b)  examined 
oblique  interactions  of  two  solitary  water  waves  of  moderate  amplitude,  using 
a  method  based  directly  on  (3)-(6)  and  which  includes  (14),  the  two-sol iton 
solution  of  (2),  as  a  special  case.  He  found  two  possible  types  of 
interactions,  depending  on  whether  or  not  two-dimensionality  dominated 
nonlinearity.  If  the  angle  between  the  two  waves  were  large  enough,  each 
segment  of  one  wave  is  affected  by  the  other  only  for  a  short  time,  and  the 
interaction  is  "weak".  These  interactions  are  qualitatively  like  linear  wave 
interactions,  as  shown  in  Figure  6a.  They  correspond  to  the  limit 
(p  -P  ) 2>> (k 2) 2  in  (14)*  Alternatively,  when  the  angle  between  the  waves 
in  small,  then  each  segment  of  one  wave  feels  the  other  wave  for  a  long  time, 
and  the  interaction  is  "strong".  Strong  interactions  are  shown  in  Figures  4, 
5,  and  6b. 
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FIGURE  6  GOES  NEXT  TO  THIS  PARAGRAPH  SOMEWHERE 


In  the  water  wave  problem,  these  strong  interactions  occur  only  if  the  angle 
between  the  two  waves  does  not  exceed  a  certain  critical  value. 
Maxworthy  (1980)  found  experimental  evidence  that  such  a  critical  angle 
exists,  and  that  the  interaction  of  two  solitary  waves  of  given  amplitude 
changes  from  "weak"  to  "strong"  as  one  moves  through  this  critical  angle. 
Johnson  (1982)  has  examined  this  critical  angle  in  more  detail. 

IV.  Periodic  Waves  in  Shallow  Water 


We  turn  now  to  the  question  of  periodic  waves  of  moderate  amplitude  in 
shallow  water,  and  to  periodic  solutions  of  (1)  and  (2).  The  derivation  of 

(I)  or  (2)  given  in  §11  remains  valid  for  periodic  initial  data,  provided  that 
the  waves  still  satisfy  assumptions  (A)-(D)  and  the  normalization  condition, 

(II) .  In  the  case  of  one-dimensional  periodic  waves,  f(r)  and  F(£)  in  (8)  are 
periodic  functions,  and  (1)  becomes  the  appropriate  evolution  equation  on 
0  <  x  <  L»  with  periodic  boundary  conditions, 

u  (x  +  L.t)  =  u  (x,t),  (15a) 


L 

J  dx  u  =  0.  (15b) 

o 

The  first  non-trivial  solution  of  the  periodic  KdV  problem  was  given  by 
Korteweg  and  de  Vries  (1895): 

u  (x»T )  *  2k2v2  cn2  [k(x-cx)  +  x0‘»  v]  +  uQ.  (16) 

Here  cn  [$;v]  is  a  Jacobian  elliptic  function  with  modulus  v  (0  <  v  <  1), 
whence  the  name  "cnoidal  wave";  also, 

c  =  6u0  -  4k2  (l-2v) , 
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kL  =  2K(v), 


u 


o 


Lhl 

K(v) 


1  +  v 


2 


> 


where  K(v),  E(v)  are  the  complete  elliptic  integrals  of  the  first  and  second 
kinds,  respectively  (cf.  Byrd  and  Friedman,  1971).  A  typical  cnoidal  wave  is 
shown  in  Figure  7.  These  waves  reduce  to  infinitesimal,  sinusoidal  waves  if 
v+0,  and  to  solitons,  (12),  if  v+1. 


FIGURE  7  NEAR  HERE 


From  one  perspective,  we  may  view  cnoidal  waves  as  the  simplest  solutions 
of  the  KdV  problem  with  periodic  boundary  conditions.  Viewed  from  a  more 
physical  perspective,  (16)  represents  a  periodic,  shallow-water  wave  according 
to 


t(x.t)  =  -|  ehu  (x »t )  +  0(e2). 

Thus,  (16)  defines  a  family  of  "typical"  nonlinear  waves  in  shallow  water, 
from  which  one  may  extract  "typical"  values  of  wave  speeds,  forces,  etc.,  for 
engineering  purposes.  This  second  viewpoint  is  common  among  ocean  engineers, 
naval  architects  and  others  responsible  for  the  design  of  large  structures  in 
relatively  shallow  water  (e.g.,  see  the  book  on  the  engineering  design  of  such 
structures  by  Sarpkaya  and  Isaacson,  1981). 

The  practical  value  of  this  approach  is  evident:  it  provides  realistic 
estimates  of  forces  from  waves  of  finite  amplitude.  Its  major  limitation  also 
is  evident:  it  is  a  one-dimensional  theory.  Water  waves  which  impinge  on  a 
vertical  wall  obliquely,  (i.e.,  waves  coming  from  almost  any  direction)  are 
excluded  from  consideration.  In  particular,  if  a  wave  impinges  on  a  wall  of 
some  structure  obliquely,  the  point  of  interaction  moves  along  the  wall  at  a 
finite  speed,  and  could  conceivably  excite  a  reasonant  frequency  of  the 
structure.  Such  wave-structure  interactions  cannot  be  predicted  using  only 
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cnoidal  waves,  because  the  one-dimensional  nature  of  cnoidal  waves  excludes 
this  possibility. 

Thus  there  is  a  practical  need  for  a  two-dimensional  generalization  of 
cnoidal  waves.  We  discuss  here  the  simplest  meaningful  generalization: 
doubly  periodic  solutions  of  the  KP  equation.  As  we  will  see,  these  solutions 
are  periodic  in  each  of  two  real,  generally  non-orthogonal  directions,  and 
except  for  degenerate  cases,  they  are  stationary  in  a  uniformly  translating 
coordinate  system. 

Without  any  calculations,  one  can  almost  guess  how  doubly  periodic  KP 
solutions  might  look.  Because  the  cnoidal  wave  is  a  periodic  generalization 
of  one  soliton,  one  should  seek  a  doubly  periodic  generalization  of  a 
two-sol iton  solution.  Thus,  Figure  8  shows  a  sketch  of  hypothetical  "weak" 
and  "strong"  interactions  of  periodic  waves,  obtained  simply  by  extending 
periodically  the  interactions  of  solitons  shown  in  Figure  6.  Note  that  with 
strong  interactions,  all  of  the  waves  are  short-crested.  As  with  Figure  6, 
one  expects  each  of  these  conjectured  interaction  figures  to  be  stationary  in 
an  appropriately  translating  coordinate  system. 

FIGURE  8  &  9  GO  ON  THE  SAME  PAGE,  NEAR  HERE. 

Preliminary  experiments  by  Hammack  (1980,  unpublished)  seem  to  confirm 
that  some  obliquely  interacting,  shallow-water  waves  form  patterns  like  those 
shown  in  Figure  8b.  The  experiments  were  performed  in  a  simple  ripple  tank 
(6.1  m.  long  x  1.13  m.  wide)  with  a  periodic  wave  maker  across  one  end  of  the 
tank.  [In  a  ripple  tank,  light  shining  through  shallow  water  from  below  is 
focussed  by  wave  crests  and  defocussed  by  troughs.  The  result  is  an 
instantaneous  picture  of  the  pattern  of  water  waves  in  which  the  wave  crests 
appear  as  bright  lines.]  In  the  experiment  in  question,  a  uniform  train  of 
one-dimensional,  shallow-water  waves  was  reflected  obliquely  by  a  uniform 
step,  placed  at  30°  to  the  incident  waves,  as  shown  in  Figure  9a.  The  water 
depth  changed  discontinuously  across  the  step  from  1.9  cm.  to  6.2  cm.  The 
incident  and  reflected  waves  constituted  two  trains  of  obliquely  interacting, 
finite-amplitude,  shallow-water  waves.  Figure  9b  shows  the  observed 
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interaction  pattern,  in  qualitative  agreement  with  that  in  Figure  8b. 
Obviously,  such  a  comparison  of  a  hypothetical  solution  of  (2)  with  a  crude 
qualitative  observation  of  water  waves  confirms  nothing,  but  it  suggests  that 
doubly  periodic  KP  solutions  might  have  practical  importance. 

Now  let  us  return  from  this  flight  of  fancy  to  the  question  of  doubly 
periodic  solutions  of  (2).  Using  methods  of  algebraic  geometry, 
Krichever  (1976)  first  observed  that  (2)  admits  quasi-periodic  solutions  of 
the  form 


u  (x.n.t)  =  2  3^  £n  ©  (Z^»  •••»  2^),  (17a) 
where  0  (Z^  ...,  ZN)  denotes  Riemann's  theta  function,  and 

Z ■  —  U.:  y  ^  Vj  t)  ^  Ui  t  ^  Z^0  ,  j=l»  •••»  N»  (17b) 
J  J  A  J  J  J  U 

(Technical  definition:  a  quasi-periodic  function  with  N  arguments  is  periodic 
in  each  of  its  arguments  separately.)  This  idea  has  been  pursued  in  a  series 
of  papers  by  Krichever  and  Novikov  (see  Krichever  and  Novikov,  1980  for  a 
review).  However,  Riemann  theta  functions  of  several  arguments  are 
mathematically  well-defined  but  poorly  understood,  and  it  has  been  difficult 
to  extract  concrete  information  from  this  work  to  date. 

Fortunately,  recent  work  by  Dubrovin  (1981)  on  compact  Riemann  surfaces 
of  low  genus  (1,  2  and  3)  has  provided  the  information  necessary  to  change 
(17)  from  an  abstract  symbol  into  a  computationally  effective  tool  for  these 
low  genera.  We  note  that  for  Riemann  theta  functions,  the  genus  of  the 
Riemann  surface  equals  the  number  of  arguments  in  the  theta  function.  Thus, 
genus  1  gives  back  the  (single-argument)  cnoidal  waves  discussed  above;  genus 
2,  corresponding  to  a  theta  function  of  two  arguments,  gives  the  intrinsically 
two-dimensional,  doubly  periodic  generalizations  of  cnoidal  waves  that  are  of 
interest  here.  We  will  not  consider  solutions  of  genus  3  in  this  paper. 
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Dubrovin's  work  enables  one  to  solve  the  following  problems: 


(I)  Parameterize  all  possible  real-valued  KP  solutions  that  can  be 
expressed  in  terms  of  a  theta  function  of  two  arguments,  as  in  (17). 

(II)  Give  an  effective  algorithm  to  compute  these  solutions. 

(Ill)  Given  the  instantaneous  values  (e.g.,  at  t=0)  of  a  KP  solution  in  the 
form  of  (17)  with  N=2,  describe  this  solution  for  all  time. 

Because  these  problems  can  be  solved,  one  can  solve  (2)  as  an  initial  value 
problem,  for  the  very  restricted  family  of  solutions  related  to  theta 
functions  of  genus  2. 

Before  describing  the  solutions  of  these  three  problems,  we  mention 
earlier  work  by  Nakamura  (1979)  and  Bryant  (1982)  on  finding  doubly  periodic 
solutions  of  KP.  Nakamura  (1979)  also  used  (17),  but  without  the  results  of 
Dubrovin  (1981)  he  was  unable  to  give  a  complete  solution  of  the  problem. 
Bryant  (1982)  gave  an  algorithm  to  compute  approximately  specific  doubly 
periodic  KP  solutions  directly  in  terms  of  Fourier  series.  No  claim  of 
completeness  was  made  in  his  work. 

(I)  Parameterize  the  KP  solutions  of  genus  2. 

Note  that  (2)  has  several  obvious  symmetries:  x_»  T>_  anc*  t-translations,  two 
coordinate  reflections  (n-*-n;  x+^X*  t-*-t),  "coordinate  rotation"  (3X+3X* 

3  -►a  +  a3  ,  3T-*-aT-6a3  -3a2a  ),  galilean  invariance  (u-*-u+e,x-*-x-66t»  n-*-n.  t+t) 

and  scaling  (u-*-A~  u,  x+*x»  n+A  n»T+A  t).  All  except  scaling  follow  from 
corresponding  symmetries  in  the  original  problem,  (3)-(6);  scaling  follows 
from  the  arbitrary  definition  of  e  in  (9).  Each  of  these  symmetries 
introduces  a  free  parameter  into  any  KP  solution.  Our  objective  here  is  to 
identify  the  parameters  that  characterize  the  KP  solutions  of  genus  2,  beyond 
these  automatic  parameters.  [According  to  this  way  of  counting,  a  single  KdV 
soliton  (12)  has  no  free  parameters,  a  cnoidal  wave  (16)  has  one  free 
parameter,  while  the  two-sol iton  solutions  of  KP  in  (14)  have  two  free 
parameters.] 


16 


Using  methods  of  algebraic  geometry,  Dubrovin  shows  that  each  KP  solution 
of  genus  2  corresponds  to  a  point  on  the  surface  of  a  compact  Riemann  surface 
of  genus  2.  (Topologically,  this  is  a  sphere  with  2  handles).  This  surface 
is  identified  uniquely  by  a  2x2  Riemann  matrix  (a  symmetric  matrix  with 
negative  definite  real  part);  i.e.,  by  three  complex  parameters.  A  fourth 
complex  parameter  identifies  the  point  on  the  Riemann  surface.  Thus  each  KP 
solution  of  genus  2  is  identified  by  four  parameters.  We  will  see  below  that 
if  these  four  parameters  all  are  real  and  satisfy  certain  inequalities,  then 
the  resulting  KP  solution  is  real. 

(II)  Calculate  KP  solutions  of  genus  2. 


where  b^j  and  Zj  are  arbitrary  complex  numbers,  Pj=0  or  1/2,  and  mj  are 
arbitrary  real  integers.  A  Riemann  theta  function  is  defined  by  a  Fourier 
series  of  the  form 

exp  |  ^  *”  *  —  *  ^  +  m  •  Z  |  ,  (19a) 


2 

I 

i  »j  =  l 

In  the  present  application  Zj  is  given  by  (17b),  and  is  pure  imaginary.  Then 
G(Z)  is  real  if  B  is  a  real  matrix,  which  we  now  assume.  £  is  negative 
definite  if 


where 


+ 

nr 


-► 

m 


bij  mi  mj 


m 


j=l 


mj  z; 


(19b) 


00  OD 

I 


V 


m  =-° 
2 


b  +  b  <0,  (20a) 

11  22 
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(20b) 


b  b  -  b  >  0. 
11  22  12 


Without  loss  of  generality  we  may  also  assume 


b  <  b  <0.  (20c) 

22  11 

We  want  (17)  to  solve  KP,  with  N=2  and  ©(Z)  defined  by  (19).  It  becomes 
necessary  to  define  certain  "theta-constants",  as  follows. 


m  =-« 
l 


®ij  cn  ■  i  S 

m  =-*  m  =-°6 
1  2 


exp  {(m+p)  .B.  (m+p)} 


(m^Pi)  (mj+pj)  exp  {(m+p)  •  | 


(21a) 


(m+p)}  (21b) 


ei  jk£ 


(mj+pj )(mk+Pk)(m£+pJl)  exp  {(m+p) 


i  •  (m+p)} 

(21c) 


For  fixed  (i,j,k,£),  each  of  these  is  a  four-component  vector,  indexed  by  the 
four  choices  of  j$. 

Define  the  4x4  real  matrix  by 

JL-  ienm  .  el2m  .  e22[«  .  eft])  .  (22) 

The  Riemann  matrix  is  said  to  be  "indecomposable"  if 

det  (D)  +  0.  (23) 

This  condition  assures  that  the  theta  function  associated  with  through  (19) 
is  honestly  genus  2,  and  has  not  degenerated  into  genus  1.  Beyond  (20)  and 
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(23),  there  are  no  further  restrictions  on  the  real  matrix  B^,  which  completely 
characterizes  the  underlying  Riemann  surface. 


Each  KP  solution  of  genus  2  corresponds  to  a  point  on  the  Riemann 
surface.  We  now  select  the  point,  by  selecting  (U  ,  U2)  in  (17b).  Because  of 
(23),  £  1  exists.  Denote  the  rows  of  this  matrix  by 


For  U  =  (U  ,  U2),  define 


a»  e  [p] 


t 


Ui  UJ  Uk  U£  Sijkt  & 


i »j  »k,£=l 


(24) 


Also  define,  for  j  >  i , 


Qij  =  X  ^  (26) 

P 

Dubrovin  shows  that  with  N=2,  (17)  defines  a  KP  solution  if  and  only  if  (U,  V, 
W)  are  related  by 

Ui  Wi  +  4  Qn  (d)  +  3  vi  =  °‘  u2  w2  +  4  Q22  (d)  +  3  v2  =  °» 

U1  w2  +  U2  W1  +  4  Q12  ^  +  6  Vi  V2  =  °-  (27) 


19 


[Note:  Dubrovin's  ft  differs  from  ours  by  a  factor  of  (-4),  because  of 
different  scalings  of  (2).]  Self-consistency  of  (27)  implies 

(U1  v2  -  u2  V2  +  |  P  (V  U2)  =  0,  (28a) 

where 

P  (U  ,  U  )  =  U2  Q  (U)  +  U2  Q  (U)  -  U  U  Q  (0).  (28b) 

If  U2^0,  the  scaling  invariance  and  coordinate  rotation  of  (2)  permit  the 
normalization. 


U  =  (s,  1),  V  =  (Vx,  0).  (29) 

Given  (29),  it  is  evident  from  (28a)  that  the  six  roots  of 

P  (s,  1)  =  0 

give  six  solutions  with  V=0,  i.e.,  genus  2  solutions  of  the  KdV  equation. 
These  solutions  correspond  to  the  six  Weierstrass  points  of  the  Riemann 
surface.  To  obtain  real-valued,  two-dimensional  solutions  of  KP,  we  need 

P  (s,  1)  >  0  for  real  s.  (30) 

Then  it  follows  from  (27)  that  (17)  defines  a  real-valued  KP  solution  of  genus 
2  if  (Z  ,  Z^)  are  defined  using 

6  -  i  (s,  1)  ,  f  «  i  (2/P(s,l)/3  ,  0)  , 

Si  =  41  (-sQ  (s,  1)  +  Q  (s,  1),  Qoo  (s,  1)).  (31) 

22  12  22 
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(32) 


There  is  one  additional  solution,  corresponding  to  U2=0.  If 


Q22  (1.0)  >  0, 


then  (17)  defines  a  real-valued  KP  solution  of  genus  2  if  (Z^  Z2)  are  defined 
using 


U  =  i  (1,0)  ,  V  =  i  (0,  2/Q22(l,0)/3) 


=  4i  (Qn  (1,0),  Q12  (1,0)) 


(33) 


This  KP  solution  is  also  a  KdV  solution  if  Q22  (1,0)  vanishes. 

The  parameter  s,  which  selects  the  point  on  the  Riemann  surface,  has  a 
simple  geometric  interpretation.  The  function  u(x,n.x)  is  defined  on  a 
fundamental  period  parallelogram,  then  repeated  periodically  in  two  directions 
(see  Figure  8).  It  is  easy  to  show  that  the  area.  A,  of  this  parallelogram 

satisfies 


(34) 


According  to  (28a),  for  fixed  scaling  of  (2),  choosing  a  solution  of  (30) 
amounts  to  choosing  the  area  of  the  period  parallelogram. 

To  summarize,  contains  three  real  parameters,  which  must  satisfy  (20) 

and  (23).  Then  s~is  a  real  parameter  that  satisfies  (30),  or  (32)  in  the 

special  case.  Every  legitimate  choice  of  these  four  real  parameters  produces 

a  KP  solution  of  genus  2,  using  (17)  with  either  (31)  or  (33).  Naturally, 

this  four-parameter  family  of  solutions  may  be  generalized  by  using  the 

synmetries  of  the  KP  equation  discussed  in  problem  (I).  All  of  the  solutions 

with  U  V  V  that  are  produced  in  this  way  are:  (i)  real-valued; 

12  2  1 

(ii)  quasi-periodic  functions  of  two  variables  (Z^,  Z2);  and  (iii)  stationary 
in  time  in  some  uniformly  translating  coordinate  system.  In  this  sense,  they 
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are  the  natural  generalizations  of  cnoidal  waves  to  two  dimensions.  We 

conjecture  that  no  other  KP  solutions  have  these  three  properties,  but  we  have 
not  proven  it. 

Figure  10  shows  some  KP  solutions  of  genus  2,  obtained  by  implementing 
this  algorithm  numerically.  All  of  the  solutions  shown  there  correspond  to 
the  same  Riemann  matrix:  bn=-1.72,  b12=1.18,  b22=-3.55.  A  phase  shift  of 

these  short-crested  waves,  due  to  their  nonlinear  interactions,  is  evident  in 
the  figures. 

FIGURE  10  GOES  HERE  OR  ABOVE. 

(Ill)  Given  the  instantaneous  values  of  a  KP  solution  of  genus  2,  describe 
this  solution  for  all  time. 

By  hypothesis,  u(x,n,T=0)  has  the  form  (17)  with  N=2.  Because  u(x,ti,0)  is 
given  pointwise,  we  may  measure  U  and  V  directly.  If  =  U£V1,  then  the 
solution  in  question  is  actually  a  KdV  solution  of  genus  2,  and  it  may  be 

described  using  the  well-established  theory  for  the  periodic  KdV  equation 
(e.g.,  Dubrovin  &  Novikov,  1974,  or  Ch.  2.3  of  Ablowitz  &  Segur,  1981).  Thus 
we  need  to  consider  only  the  case  in  which 

V2  *  u2Vi  •  .  (35) 

so  that  the  solution  is  non-trivially  periodic  in  two  spatial  directions. 

Every  KP  solution  of  genus  2  that  satisfies  (35)  is  stationary  in  a 
uniformly  translating  coordinate  system,  whose  velocity  is  given  by  W  in 

(17b).  Once  W  is  known,  then  the  initial  data  in  a  period  parallelogram  plus 
W  determine  the  solution  fcr  all  time.  It  is  not  necessary  to  reconstruct  j5, 
the  Riemann  matrix,  for  KP  solutions  of  genus  2  that  satisfy  (35). 

Algebraic  equations  for  Cl  may  be  obtained  in  a  variety  of  ways.  The 
method  presented  here  is  valid  either  if  U2=0  or  if  (U x /U2 )  is  rational.  This 
restriction  is  always  satisfied  in  applications,  where  (Uj/U2)  is  measured 
only  to  a  finite  accuracy.  In  problem  (II),  it  amounts  to  requiring  that  the 
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parameter  s  be  rational. 


If  (Ul/U2)  is  rational  (or  if  U2=0),  then  u(x»n,0)  is  a  strictly  periodic 
function  of  x»  holding  (n *t )  fixed.  Denote  this  x-Period  by  L*  Because 
uCxt’i.O)  has  the  form  (17),  it  follows  that  for  all  (xQ»n). 


0. 


(36) 


Defi ne 


^(x.nJx^)  =  ~  i  J 


d5  •  5u(x+x  +£.’■)). 
o 


(37) 


4>  is  the  unique  anti-derivative  of  u  with  the  same  periodicity  as  u,  and  with 
zero  mean  in  x>  1.e.»  4>  also  satisfies  (36),  and  3^4>=u.  The  corresponding 
anti -deri vati ve  of  <f>  may  be  defined  in  a  similar  way. 


One  equation  for  ft  may  be  obtained  by  multiplying  (2)  by  u,  and 
integrated  over  a  period  parallelogram  (with  area  A).  The  result  is 


dxdn  [u^u 


T 


+  6uu 


u2  +  3u2]  =  0  . 
XX 


(38) 


By  hypothesis,  u(x»n,T)  =  f(Zj,  Z2),  where  Zj  is  given  by  (17b).  It  follows 
that 

(U1V2  ‘  Vl>  UT  ’  (“lV2  ‘  W  UX  +  (U1W2  '  W  “n  •  (39) 

so  (38)  becomes 
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// 


dxdn  [V2u^  -  U2  uxun]  W 


A 


A 


ff 

J  J  dxdn  [6uu 
A 


X  XX  1 


2  2  2 

-  u  +  3u  ]  =  0  . 


(40) 


To  obtain  a  second  equation,  multiply  (2)  by  that  anti-derivative  of  $ 
(i.e.,  the  second  integral  of  u)  with  the  same  periodicity  as  u  and  $.  The 
result  after  integrating  over  the  same  period  parallelogram  is: 


A 


A 


(l)  V  -  U  V  ) 
v  1  2  2 


(41) 


A 


Equations  (40)  and  (41)  are  two  linear  algebraic  equations  for  (W  ,W  ). 
Wherever  they  are  linearly  independent,  their  common  solution  defines  (W  ,W2), 
which  completes  the  mathematical  specification  of  the  KP  solutions  of  genus  2. 
It  is  easy  to  show  that  (40)  and  (41)  are  linearly  independent  for  small 
enough  wave  amplitudes,  but  we  have  not  yet  established  this  property  in 
general.  In  the  cases  we  have  tested  numerically,  (40)  and  (41)  are 
independent,  and  their  common  solution  agrees  with  that  in  (31). 

More  generally,  given  any  initial  data,  f(x,n)»  such  that 


L 


o 


(42) 
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one  may  seek  the  KP  solution  of  genus  2  that  best  approximates  f(x.Ti)  at 
t  =  0,  along  with  an  appropriate  measure  of  goodness  of  fit.  It  is  evident 
that  such  approximation  procedures  will  be  necessary  if  these  KP  solutions  of 
genus  2  are  to  become  a  practical  tool  in  physical  problems.  However,  we 
leave  this  and  other  questions  of  physical  implications  of  these  solutions  for 
a  future  paper. 

We  are  grateful  to  T.  Toedtemeier  and  J.  L.  Hammack  for  their  permission 
to  use  Figures  5  and  9,  respectively.  We  acknowledge  with  pleasure  several 
helpful  conversations  with  R.  I.  Sykes  concerning  the  algorithms  in  §  4.  This 
work  was  partially  supported  by  the  Office  of  Naval  Research,  and  by  NSF  Grant 
#MCS-8108814  (A01 ). 
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Figure  1 
Figure  2 

F  i  gu  re  3 

Figure  4, 
Figure  5. 
Figure  6. 

Figure  7. 
Figure  8. 

Figure  9. 

Figure  10. 


Figure  Captions 

Physical  configuration,  showing  notation  for  Eq.'s  (3)-(6). 
(1974)!1C  drawing  0f  the  wave  9enerator  used  by  Hammack  and  Segur 


Evolution  of  surface  waves,  measured  at  four  locations  downstream 

?!  *ave  maker*  The  front  of  each  wave  is  to  the  left  in  the 
coordinate  system  used,  which  translates  with  speed  /gh.  _ 

*ave  Profiles;  .  soliton  profiles  computed  usinj 

(13)  and  the  measured  peak  amplitude  of  each  wave.  ^  9 


Two  soliton  solution 
4  in  (14). 


Pl  =  -p2 


of  the  KP  equation,  with 


kr 


k2  - 


i. 


Oblique  interaction  of  two  waves  in  shallow  water  observed  off  a 
beach  in  Oregon.  (Photograph  courtesy  of  T.  Toedtemeier. ) 

patterns*  0°/  interactions  of  Kp  solitons,  showing  the 

patterns  of  the  wave  crests.  Arrows  indicate  direction  of 

anol’e9hPi!0n‘  A*  Weak  interaction.  corresponding  to  a  large 
angle  between  the  waves,  or  to  small  amplitude  waves*  (b)  strono 

ImpmSde^avIs!'  8  S"a11  a"9,e  bet“ee"  the  wa,es-  “r  t0  ,ar9er 


One  period  of  cnoidal  wave 
v2  =  1/2  in  (16). 


solution  of  the  KdV  equation,  with 


Hypothetical  pattern  of  wave  crests  for 
(a)  Weak  interaction,  as  in  linear 
interaction.  The  dashed  lines  indicate  a 


periodic  KP  solutions. 

theory;  (b)  strong 
period  parallelogram. 


(a)  Experimental  apparatus  to 
trains  of  finite  amplitude  waves 
configuration  of  wave  crests. 
Hammack. ) 


create  oblique  interaction  of 
in  shallow  water,  (b)  Observed 
(Photograph  courtesy  of  J.L. 


Sample  of  KP 
points  on  the 
(b)  s  =  0.35, 


solutions  of  genus  2, 
same  Riemann  surface, 
(c)  s  =  -0.70. 


corresponding  to  different 
(a)  s  =  Uj/Uj  =  0.0  in  (29); 
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